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Abstract 

We study shell models that conserve the analogues of energy and enstrophy, 
hence designed to mimic fluid turbulence in 2D. The main result is that the 
observed state is well described as a formal statistical equilibrium, closely 
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analogous to the approach to two-dimensional ideal hydrodynamics of On- 
sager, Hopf and Lee. In the presence of forcing and dissipation we observe 
a forward flux of enstrophy and a backward flux of energy. These fluxes can 
be understood as mean diffusive drifts from a source to two sinks in a system 
which is close to local equilibrium with Lagrange multipliers ("shell tempera- 
tures") changing slowly with scale. The dimensional predictions on the power 
spectra from a supposed forward cascade of enstrophy, and from one branch of 
the formal statistical equilibrium, coincide in these shell models at difference 
to the corresponding predictions for the Navier-Stokes and Euler equations in 
2D. This coincidence have previously led to the mistaken conclusion that shell 
models exhibit a forward cascade of enstrophy. We also study the dynamical 

properties of the models and the growth of perturbations. 
47.25.C 
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I. INTRODUCTION 



The idea of shell models of turbulence is to replace the fluctuations of a turbulent 
field in an octave of wave numbers 2 n < \k\ < 2 n+1 by one or a few representative variables. 
The range of wave numbers is called a shell, and the variables are called shell variables. The 
dynamics of the shell variables should be chosen so as to preserve as many as possible of the 
qualitative features of the full equations. In this way, it is possible to describe the cascade 
processes of inertially conserved quantities by a chaotic dynamical system with a limited 
number of degrees of freedom. Predictions from dimensional analysis can then be tested 
with far greater accuracy than in full simulations of the Navier-Stokes equations. 

This approach has been particularly successful for 3D turbulence in the models intro- 
duced by Gledzer f8[ and Yamada & Ohkitani || (see also || or |§), which, following the 
recent literature, we will refer to as the GOY models. Nevertheless it has recently become 
clear that the situation is more complicated. Particular variants of the GOY models have 



in fact stable fixed points, corresponding to the Kolmogorov scaling law in 3D [ 10 1 . The 



dimensional predictions for the energy spectra come out correctly in these models ||10|| , but 
there are no intermittency corrections, and the whole phenomenology of a turbulent state 
(e.g. sensitivity to initial conditions, positive Lyapunov exponents) is absent. 

We will here report on another instance where shell models give significantly different 
results than expected at first: shell models designed to model fluid turbulence in 2D. We 
will show that these models are perfectly well described by a formal statistical mechanics, 
closely similar to the approach of Onsager |Tl| , Hopf [12|] and Lee [0] for inviscid 2D 
hydrodynamics (see also the review [|14[)- If energy and enstrophy are pumped into the 
system by an external force and removed by viscous terms, we do observe a net mean flux 
of both energy and enstrophy from the force to the viscous sinks. These fluxes are shown to 
be analogous to mean drifts proportional to the gradients of the conjugate quantities, as in 
a system which is locally, but not globally, in thermodynamic equilibrium. There is no need 
to invoke a cascade, neither in the inverse direction (of energy) nor in the forward direction 
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(of enstrophy). We show that in these shell models the dimensional predictions on the 
energy spectrum from a supposed forward cascade of enstrophy coincide with equipartition 
of shell enstrophies, which is one of the branches of the statistical equilibrium. In previous 
investigations [|i~5|-|T7|j the simpler explanation of the observed spectrum in terms of a weakly 
perturbed equlibria was overlooked. 

The paper is organised as follows: in section |T| we review the standard phenomenology of 
2D turbulence and the minimal constraints which have to be put on shell models to take into 
account the difference to 3D. We present the shell models tested in the simulations which 



are all extensions and variations of the 2D GOY model studied by several authors 18.15-17 



In section III we recall the statistical equilibrium approach, and apply it to determine the 



scaling laws and the probability distribution functions of the shell models. In section [TV 
we present numerical results on a steady forced state with enstrophy and energy output, 
and compare them with the predictions of section |HI| . In section [V| we present the results 



on dynamical properties, which are quite different from the results for 3D shell models, 
but consistent with what is expected from a system close to equilibrium. In section 6 we 
summarize our results. 

II. 2D TURBULENCE AND SHELL MODELS 

The main particularity of 2D Euler equations is that the vorticity is a material invariant, 
i.e. conserved along stream lines. This leads to an infinite number of additional integrals 
of motion, since the global integral of any functional of vorticity is inertially preserved. In 
practice one habitually considers the lowest polynomial positive invariant, the enstrophy. 
One of the main reasons for this is that the enstrophy is also an invariant of the truncated 



Euler equations in Fourier space |TJj]. Although energy and enstrophy are both conserved 
globally, their support in wave vector space may change over time. The changes must be in 
opposite directions, such that if enstrophy is transported towards small scales (large wave 
numbers) energy is transported to large scales (small wave numbers). 
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The Navier-Stokes equations describe in addition to the inertial motion the action of 
external forcing and molecular viscosity. The first major difference to the situation in 3D is 
that in 2D molecular viscosity cannot remove energy efficiently, because enstrophy contains 
two more derivatives, is dissipated first, and energy is then constrained to flow to large 
scales. This implies that in numerical investigations an infra-red viscous term v\— A) -7 is 
necessary to remove energy at large scales: 

d t u+ {u-V)u = -V P + uAu - v \-A)- J u + F (f) 
V-n = V-F = 

The exponent 7 is rather arbitrary. For instance, 7 in (1) equal to zero corresponds to linear 
friction, which describes the physical situation of a fluid flow in thin films with viscous drag 
in boundary layers. If one is mainly interested in the backward energy cascade ranges it is 
better to take a sharper artificial infra-red viscosity, for which the action is more concentrated 
in a small range of wave number. This procedure is standard in numerical investigations of 
the 2D Navier-Stokes equations (see e.g. P^||TP|JT^] ). 

We will from now on assume that the force acts only on scales with typical wave number 
kf, and that the molecular and infra-red viscosities act on scales respectively much smaller 
and much larger. The spectrum can then be divided into five ranges: the energy dissipation 
range, where the infra-red viscosity acts effectively; the inverse inertial transport range, 
where energy flows towards large scales; the scale of the forcing; the forward inertial transport 
range, where enstrophy flows towards small scales; and the enstrophy dissipation range, 
where molecular viscosity acts effectively. 

The first attempts to predict the behaviour of hydrodynamics in 2D used a formal statis- 
tical mechanics analogy for the Euler equations. We review this approach in section PTT| . A 
later approach due to Kraichnan and Batchelor proceeds along the lines of the Kolmogorov 



cascade picture for 3D ||20|| . For the forward transport range, a cascade of enstrophy implies 
by dimensional arguments that the energy spectrum should be E(k) ~ k~ 3 . Similarly, for 
the inverse transport range, an inverse cascade of energy gives E(k) ~ fc-s Hm|l] 



In 3D there is strong experimental support of the cascade picture [TJI| , albeit with inter 



mittency corrections |2^J2^] . The evidence for cascade processes in 2D is not as clear. In the 
forward range, a energy spectrum steeper than k~ 3 implies that the most important inter- 
actions are not local in wave number space. A scaling k~ 3 is marginal, and the dimensional 
predictions based on a cascade process with local interactions in fc-space are questionable. In 
fact, most numerical investigations of the 2D Navier-Stokes equations report spectra which 
are significantly steeper than k~ 3 [^6,^7]. However, all numerical investigations of the full 



Navier-Stokes equations suffer from a limited range of wave numbers for the inertial ranges. 
The k~ 3 energy spectrum can therefore not conclusively be ruled out, while the prediction 
of the statistical equilibrium approach, which is k~ l , is clearly excluded by the numerical 
observations. 

The true state of the 2D Navier-Stokes equations in the inverse cascade range is un- 
clear. Most simulation have reported a k~ 3 spectrum in agreement with the dimensional 
predictions fl8| . |2"8|j . The latest investigation over much longer time-scales instead reports a 



k 3 spectrum extending to scales about one order of magnitude larger than the forcing, and 
then a spectrum clearly much flatter than k~i |T9|]. In the inverse range the prediction from 



statistical equilibrium, which is E(k) ~ k, does therefore not seem to be completely ruled 
out. 

The shell models we consider should be taken as the simplest systems that share the 
property of having two inertially conserved quadratic quantities. They are systems of coupled 
ordinary equations which may in general be written as: 

jW n = J2 XnmlWnWt - K n W n + F n (2) 
m,l 

where X indicates the nonlinear interactions, K the viscous forces and F an external force. 
In (Q) we have suppressed an index of shell variables within one shell. In the extreme case 
considered in the literature J7| this can be a moderately large number, of the order of 10 2 . 
In the cases we consider, we will have one or a few complex shell variables per shell, with 
nonlinear terms that are compactly expressed using the complex conjugates W*. 
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We use the convention that the shell variable W n is like a shell vorticity, so (§) is a model 
of the vorticity equation. The quantities that should be inertially conserved are then: 

"enstrophy" ; Q = £ | W n | 2 (3) 

n 

and "energy"; E = EIt^| 2 > (4) 

n ™n 

where k n is the typical wave number of shell n. We will only consider the case where the 
fc n 's are spaced in octaves, i.e. k n = k Q 2 n . 

The dimensional predictions are based on the assumption of cascade processes and from 
the observation that the fluxes of energy U n (E) and of enstrophy IT n (f2) through the n th 
shell are given by certain third order correlation functions 0. For the 2D GOY model (see 
equation (|T2|)) they read in detail 

U n (Q) = (Re(4W n W n ^W n ^ - W n + x W n W n _ x )) (5) 
U n (E) = ^(Re(W„W n _ 1 W n _ 2 - W n+1 W n W n ^)) (6) 

If there is a net forward transport of enstrophy from the force to the enstrophy dissipation, 
the correlation functions (|5|) (and the analogous expressions for the other models) must be 
constant in the forward transport range. That gives the dimensional estimate 

(\W n \) ~ 4k° n , (7) 

where is the mean dissipation of enstrophy per unit time. Similarly in the inverse transport 
range one gets the estimate 

(|Wg 2 > ~ elkl, (8) 
where e is the mean dissipation of energy per unit time. The shell energy estimates are 

_ 2 

respectively E n ~ k~ 2 and E n ~ kn 3 , which is the same as the cascade estimates E(k) ~ k~ 3 
and E(k) ~ k~% for the full Navier-Stokes equations integrated over one octave in wave 
number space. 
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A. The GOY models in 2D 



The GOY models are usually written in terms of complex shell velocities as 

jU n = ik n {a n U: +l U: +2 + + ^U:_,U:_ 2 ) - uk 2 n U n + /„ (9) 

We can transform (Q) to (Q) by absorbing a factor k n in the shell variables. The particular 
form of the right-hand side of (|9J) ensures that the phase-space volume 

dV = dUi A dU* A ... A dU n A dU* (10) 

is inertially conserved. The parameters values for which (f|) conserve energy are given up to 
an overall factor by 

a n = 1 b n = -e c n = -(l-e) (11) 

The bifurcations and transitions to chaos in (|]) when changing e have recently been inves- 
tigated in JK|. If e is greater than 1 there is a second conserved quantity; J2kn a \U n \ 2 , the 



exponent a changing with e 



Enstrophy conservation, i.e. a = 1, is realized at the particular value e equal to 1.25. 
The 2D GOY equations written for shell vorticity variables hence read 

jW n = (w: +1 w: +2 - sC-iC+i + ^UWU) - <w n - v'K 2 m n + F n (12) 



For convenience we have in (|12D rescaled the W n 's to absorb an overall factor |, and added 
an artificial infra-red viscosity. 



In [0 a class of models were proposed that contain GOY-like interactions between triads 
of shells further apart. These models were further studied in the range of forward transport 
in [|17]]. For interactions decaying sufficiently fast with distance between interacting shells 
they behave quite similarly to the 2D GOY model. 
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B. The coupled GOY model 



We motivate this model by the observation that the GOY models make no difference 
between velocity and vorticity, which only differ by a scale factor. It could therefore be 
hoped that a shell model which preserves some trace of the vector structure of the velocity 
field will be in a qualitative sense closer to the Navier-Stokes equations. 

In 2D we can write the vorticity equation in Fourier space as 

j t m= £ ^^u(P)u(P')-u\k\"u(k) + F(k) (13) 

k+k'+k"=0 l fc I 

An obvious consequence of ( |T3"D is that two parallel vectors do not drive a third. We can 
model (|i~3"D by taking few shell variables, W n j, j = 0, 1, . . ., per shell, associate to each a 
wave vector k n j with length k n and direction e.j, and look for equations in the form 

(14) 

-vk 2 n w nJ - v'K 2 m n , 3 + F nd 

with some interaction coefficients C chosen to preserve energy and enstrophy. 

A simple implementation of (13) is to take wave vectors in a hexagonal pattern, 



ij = (cos(27r/'/3), sin(27r//3)) for j = 0,1,2. One possible set of interactions are then 
the same as in the 2D GOY model, but going only between triples of shell variables having 
different directions: 

iwnj = m +1>j+1 w* +2>j+2 - 5w: +w w:_ 1J+2 + m;_ 2tj+1 w;_ 1>j+2 ) + 

m +w w: +1J+2 - hw*_ ld+1 w: +ltj+2 + w:_ 1J+l w:„ 2d+2 ) (15) 

-uk$W n j - v'k~ 27 W n j + F nij 



In ( P~5| ) we have used the notation that W n j + z = W n j. We will call ( [ToD the coupled GOY 
model. It has some moderate numerical advantages over the original GOY model. Due to 
more variables per shell and more couplings, the coupled GOY model equilibrates faster. In 
the original GOY models with a steady force one observes that the shell energy spectrum 



in the inertial range exhibits oscillations superimposed on a mean power-law ||. We find 
no trace of such oscillations in (^5|), and conclude that three variables per shell is enough to 
sufficiently randomize the system and remove these undesired oscillations. 
We have used forcing terms acting on shell of the forms: 

F "-> = ^rki- w "-> (16) 

3=0,1,2 

Foj = VkiF 1 ^- (17) 

3=0,1,2 

The advantage of ([16]) and ([17D is that the fluxes of energy and enstrophy from the force are 
fixed to be proportional (in the case (|T6| ) equal) to t]. For most of our simulations we have 
used artificial viscous terms like —v'k~ 2 W n , i.e. 7 = 1. 

III. THE STATISTICAL MECHANICS APPROACH 



We here review the approach from equilibrium statistical mechanics of Onsager |TT 
Hopf Jl2| and Lee |13 ] . For a later review see jbf . 



We consider the equations of two-dimensional hydrodynamics in Fourier space, with cut- 
offs k min and fc max . The lower cut-off can be realised physically by enclosing the system 
in a finite container. We need also the upper cut-off to avoid the ultra-violet catastrophe 
of classical continuous fields. As the system leaves inertially invariant the two quadratic 
forms energy and enstrophy, the canonical ensemble distribution function for the Fourier 
components of vorticity is 

P{u>(k)) oc exp[ - (J3i \u(k)\ 2 + (3 2 \k\- 2 \u(k)\ 2 )/2] (18) 

In ( [IB] ) we have two Lagrange multipliers, f3i and 02, or, alternatively, an "enstrophy tem- 
perature" Ti, equal to l//5i, and an "energy temperature" T 2 , equal to l//?2- The second 
order moment of ( |IS| ) is 

<^>i 2 > = irf^s (19) 
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which, depending on \k\ 2 and the temperatures, separates into two branches 

\t 2 \k\ 2 if \k\ 2 < ^ 

(\m\ 2 ) ~ 'J Ta (20) 

[ Ti if \k\ 2 » | 

The energy spectra in these two ranges are E(k) ~ k and E(k) ~ respectively. As 
discussed above at least the second prediction is clearly in disagreement with all numerical 
investigations of the 2D Navier-Stokes equations |18| , p8| , p!9| , p6 . 



The argument for (|T8|) translates immediately to the shell models considered in section ||. 
We then have the prediction that the shell variables are uncorrelated Gaussians with widths 

[ T 2 k\ if k\ < ^ 
(\W n \ 2 ) ~ (21) 
[ T x if A* » | 

The two branches in (|2T| ) simply correspond to shell energy equipartition and shell enstrophy 
equipartition. The shell energy prediction in the shell enstrophy equipartition branch is 
E n ~ k~ 2 . We thus have that in 2D shell models, but not in 2D ideal hydrodynamics, the 
energy spectrum of one branch of a formal statistical equilibrium and of an assumed forward 
cascade of enstrophy exhibits the same scaling law. 

We continue with some elementary considerations from non-equilibrium statistical me- 
chanics Suppose that a physical system is not in thermal equlibrium, but the distribu- 
tion functions at any point are close to the equilibrium distribution, with temperatures that 
change slowly in space. A heat flux is then set up, which acts to restore global equilibrium. 
In other words, the induced flux is 

n Q = -kVT (22) 

where the proportionality constant k is the heat conduction coefficient, generally dependent 
on temperature. 

We now translate this picture to the shell models. For defmitiveness we consider the range 
of forward transport of enstrophy. The local "shell temperatures" of this non-equilibrium 
system are not defined per se, but may be specified by the expectation values 
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Ti(n) = (\W n \ 2 ) (23) 

where the average is taken with respect to the the stationary distribution of the W n 's. In 
first approximation, it is assumed that this joint distribution factorizes into uncorrelated 
Gaussians with widths T\(n). Then, if the shell temperatures of two nearest neighbour 
shells differ, we expect a flux of enstrophy through the n th shell. 

U n (Q) = -KtiT^n + 1) - Ti{n)) (24) 

where we can call the transport coefficient K\ the "shell enstrophy conductivity constant". 
In a steady state of forward flux of enstrophy, one thus has 

T 1 (n+l)=T 1 (n)-^ (25) 

where e w is the mean dissipation of enstrophy per unit time. 

The temperature dependence of K\ can be estimated as follows. The nonlinear inertial 
term in the shell model equations plays the role of the collision operator in kinetic theory. 
The transport coefficient is determined by the linearized collision operator. In the shell 
models the nonlinear terms are quadratic, and the linearization hence linear in the shell 
amplitudes. We therefore expect 



K\ ~ 0y7i (^) 

with some proportionality constant a, as in the kinetic theory of gases. 

When all the temperature increments are supposed small we have the prediction 

T i( n ) ~ [( T i( n f))^ ~ 7^(« - n/)]s (27) 

where rif is the shell of the forcing. At the dissipative end T\{n) eventually becomes much 
smaller than unity. We may therefore rewrite (p7| ) as 



TWn) 



/ \ 

Za 
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We see that by coincidence the dependance on is also as predicted by the cascade picture. 
The spectrum however depends logarithmically on v through n^ ss . This contradicts the 
assumption of a cascade process independent of viscosity. 

Finally, over small increments the transport coefficient are slowly varying, and we have 
the shell model equivalent of Fourier's law of heat conduction: 

Ti(n) ~ Ti(n ) - — (n - n ) (n - n ) < (n diss - n f ) (29) 



IV. NUMERICAL RESULTS ON THE STEADY STATE 

In this section we present numerical results on the coupled GOY model (1151) . The first 
results, figures [1] to |6| are from one run of length 3 • 10 5 in time. Statistics was taken with 
frequency 1 Hz, about once per shell turnover time in the range of forward transport. By 
convention the force acted on shell zero. The viscosity was v = 10~ 24 . Practically all 
dissipation of enstrophy occurred in shells 34 to 42. The artificial viscous term was of the 
type —v'k~ 2 W n with v' = 10~ 15 . The energy output occurred in shells —21 to —14. We used 
the forcing of type (|i~7|) with rj equal to one. The observed mean fluxes of enstrophy and 
energy were approximately 0.075. The integration method was the slaved leap-frog scheme 
of !| with time-step 0.003. 

Let us first look at figure [l] which shows the second moment of shell vorticity vs. shell 
number. The dominant overall feature is one branch at positive n's which is nearly flat, i.e. 
implies (l^nl 2 ) ~ and one branch at negative n's which is closely fitted by (IW^I 2 ) ~ /c 2 . 



Both results are in agreement with the predictions of a statistical equilibrium as in section pi. 
Let us add here that according to figure [1] the characteristic times in the range of forward 
transport of enstrophy are about constant. This run therefore represents about 10 5 turnover 
times in all the forward range, and there are indeed only small fluctuations. In contrast, 
in the inverse range the characteristic times are proportional to k" 1 , which means that we 
have about one turnover time in the energy dissipation range. Looking closer on figure [IJ 
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we clearly see that the spectrum is in the range of forward transport is not quite flat, and 
that the "fc°-range" seems to extend about five shells to the left of the force. The second 
phenomenon is difficult to explain in the framework of the cascade picture. It is not a 
problem in the statistical equilibrium picture if we assume that the shell temperatures of 
enstrophy and energy are such that the bend in the curve occurs somewhat to the left of the 



force. The first phenomenon was observed in ||17|| , and there attributed to a non power-law 



correction on top of a power-law with an exponent close to zero. Systematic investigations 



on the Yamada-Ohkitani models [|3l] gave only a pre-asymptotic correction, and thus an 
asymptotic power-law k^. Figure H shows the second moments of the shell variables in real 
scale, instead of logarithmic scale as in figure [1]. A dominant linear trend is quite clear. The 
deviations from a pure power-law in figure [I] are hence as predicted by (^9]) from the picture 
of weakly perturbed equilibria. 

Figure |3] shows the numerically observed distributions in the forward range, compared 
with a Gaussian. The fit is clearly good. Figure ^ compares the numerically observed 
kurtosis of the shell enstrophies to the estimate that the three variables in one shell are 
uncorrelated Gaussians with the same variance. Again the fit is quite good. This analysis in 
the inverse range leads to similar scenario. Taken together these results imply that to good 
approximation there are no observable correlations between the shell variables within one 
shell. We have also checked directly that this is true. In addition there are no measurable 
correlations between shell variables in two nearest neighbour shells. 

Figure [5] shows the mean enstrophy flux and the mean squared enstrophy flux. The curve 
of mean enstrophy flux should be flat, and is so to a much better approximation than the 
second moment in figure |l|. 

Figure |5] shows that the standard deviation of the flux is about 50 times higher than the 
mean value in the beginning of the range of forward transport, and then decreases, but stays 
above the mean flux all up to the dissipative range. The mean flux shows a clear downwards 
jump from shell one to shell zero, and then a fast decrease for negative n's. In this range 
the eventual mean transport can only vanish, and the non-zero values represent fluctuating 
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incremented averages taken over finite time. 

Figure |6| shows the mean energy flux and the mean squared energy flux. Towards the 
energy dissipation range the mean flux is to good approximation constant, and agrees with 
the mean energy infused to the system per unit of time. At the top value around five shells 
to the left of the forcing, the standard deviation of the energy flux is about 10, 000 times 
larger than the mean value. From the value of the variance of the flux, about 10 6 , and the 
number of measurements, 3 x 10 5 , one sees that rise of the curve of the energy flux below the 
top value of the flux squared curve is a statistical measurement error. We have checked that 
this explanation is correct by comparing with other runs with higher artificial viscosity for 
which the range of inverse energy transport is shorter, and the times to reach an asymptotic 
state at the energy dissipation scale is shorter. 

The second result are direct tests of equation (p6|), from different runs with different 
viscosities and different forward enstrophy flux. Figure [7| shows that K\ scales to good 
approximation as y^T as expected by our theoretical arguments. 



V. DYNAMICAL BEHAVIOUR 

In this section we briefly discuss the "instability properties" of the 2D shell models flT5D 
by studying the behavior of the Lyapunov exponent A as function of the Reynolds number 
Re, the fluctuations of the effective Lyapunov exponent j T , the tangent vector and the 
spatial spreading of a small perturbation initially localized on a given shell. 

Given a dynamical system described by a set of differential equations 

% = F[x(t)} (30) 

the response of the system to a perturbation Sx(t + r) of its state at time t after a delay r 
is measured by the error growth rate 

ft( T ) J*+J>Ufe<' + J>l (31) 
\z(t)\ \Sx(t)\ 

where z is the tangent vector obeying the evolution equation 
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dzi _ ^ dFj[x(t)] 

IF \ ~~iJ7, {6Z) 

By definition the largest Lyapunov exponent is 

A = Km - (In i2 t (r)) (33) 

where the angular brackets denote a time average along the trajectory. The Oseledec theorem 



32[ ensures that if the average is removed, then for almost all initial conditions one obtains 
the same value for A. 

The exponent A gives a global characterization of the "instability" of the trajectory. 
Local informations can be obtained from the effective Lyapunov exponent j T (t) defined as 

7r(f) = - In Rt{r) (34) 
r 



and from its fluctuations p3 | 



The scenario which emerges is completely different from that observed for the 3D GOY 
model studied in ||. The main results for the 3D model are that: 

1. the Lyapunov exponent A increases with the Reynolds number: 

A ~ Re a with a ~ 0.46 (35) 



in good agreement with the prediction of the multifractal generalization [34] of a Ruelle 



argument |35| 



2. the effective Lyapunov exponent 7 r exhibits strong fluctuations, at increasing Re. Its 
variance fi scales like 

H = hm t({j t - A) 2 ) ~ Re p with (3 ~ 0.8 (36) 

so that the ratio /i/A, which gives a quantitative measure of the intermittency level 
diverges with the Reynolds number. 
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3. the tangent vector during the intermittent bursts (of energy and chaoticity) is strongly 
localized on the shells corresponding to the Kolmogorov scale, while in the laminar 
phase spreads over all the inertial range. In addition there is a strong correlation 
between the intermittent bursts of energy dissipation, large fluctuations of the effective 
Lyapunov exponent and the localization of the tangent vector on the Kolmogorov 
length. 

4. there is a backward cascade for the propagation of a perturbation from small scales to 
large scales in qualitative agreement with the phenomenological scenario proposed by 



The 2D shell models present a completely different behaviour. The numerical study of 
the models introduced in this paper reveals the following scenario: 

1. the Lyapunov exponent depends very weakly on the Reynolds number: 



see figure [|. It is worth noting that even in the framework of the Kraichnan-Batchelor 
theory, one has the same prediction for 2D Navier-Stokes equations, by assuming 
that the Lyapunov exponent is proportional to the smallest characteristic time of the 
system. 

2. the effective Lyapunov exponent has small fluctuations and /i < A. 

3. the tangent vector is concentrated at the scale of the forcing. 

4. a small perturbation initially concentrated about a given shell does not propagate by 
an inverse cascade mechanism, but diffuses through the shells. In Fig. ^ we report the 
spreding of a perturbation initially located on a small number of shells. The quantity 
shown in figure is defined as 



Lorenz 36 



A ~ In Re 



(37) 



1 V2 




(38) 



i=o,i,2 
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where at the initial time W n ^ differs from W n j of a quantity 5 on shells n = n\, . . . , n%. 
In the case reported in Fig. |9| we used n\ = —6, n 2 = —5 and 5 = 10~ 10 . However the 
curves are rather insensitive to the initial location of the perturbation. 

The above results provide a clear evidence the "instability" behavior in the 2D shell 
models is very different from those of 3D shell models. Let us remark that the dynami- 
cal features observed for equation ([15]) are rather close a scenario expected in equlibrium 
statistical mechanics. 

VI. DISCUSSION 

We have presented numerical results on a wide class of shell models of 2D turbulence 
which can be simply and coherently explained by a formal non-equilibrium statistical me- 
chanics, close to local equilibrium. The presence of an inverse cascade of energy can be ruled 
out already by analyzing the second moments of the shell variables. In the forward inertial 
range a direct cascade of shell enstrophy would predict the correct scaling behaviour of the 
second moments with shell wave number and mean dissipation of enstrophy. However, we 
also find that the spectrum rises as the viscosity decreases. This rules out a cascade process 
which assumes a state independent of viscosity in the low viscosity limit. 

The weakly perturbed equilibria predicts Gaussian probability distribution functions of 
the shell variables, and the pre-asymptotic corrections to the power-law in the forward range. 
Let us remark that these two observations are difficult to reconcile in a cascade picture where 
Gaussian probability distribution functions are possible if there is only one scaling exponent, 
"unifractality" , ||37|| . However, in absence of multifractality, it is hard to imagine mechanisms 
which are able to produce corrections to scaling such as pseudo-algebraic power laws, the 
so-called multiscaling |38fl . 

Finally, the mean values of the fluxes of both energy and enstrophy are always much 
smaller than the standard deviations. The fluxes are thus always small corrections superim- 
posed on an mean randomly fluctuating state. 
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For the full hydrodynamic equations in 2D a state of equilibrium is ruled out in the range 
of forward transport. We therefore conclude that the models studied in this paper have very 
little to do with 2D turbulence. 

A natural question at this point is: why have shell models for 3D turbulence given 
reasonable results, while the 2D models have not? A qualitative answer goes as follows: 
the statistical equilibria picture should be relevant if the time-scales of relaxation to local 
equilibrium are faster than the timescales of transport of the conserved quantities to the 
viscous sinks. The time-scale for relaxation to local equilibrium can be estimated as the local 

-2 

shell turnover times, which in 3D shell models decrease with shell number (as t n ~ kn at 
the Kolmogorov fixed point). The time scale for transport to the energy dissipation range 
from a shell in the forward transport range can then be estimated as the geometric sum of 
the local turn-over times up to the energy dissipation shell. Hence both time scales are of 
the same order, and it is unlikely that 3D shell models display statistical equilibrium. And 
indeed a cascade scenario is observed. 

On the contrary, in 2D shell models the local turnover times in the forward range are 
constant, and the time to transport enstrophy to the dissipation range is proportional to 
how far away in shells that is. Therefore local statistical equilibrium has a chance to develop. 

The question then arises why this argument does not work for the 2D Navier-Stokes 
equations, where, in the Batchelor-Kraichnan cascade picture, the timescales in the forward 
range are also constant. One possible explanation is that one has to take into account the 
full nonlinearities of the Navier-Stokes equations with the non-local transfer of energy and 
enstrophy. A more interesting possibility is that the reason is that there are more states 
at high wave-numbers, and that the system gains entropy locally by transporting in that 
direction. It would then make sense to study simplified models but where the number of 



states increases with shell number. Numerical experiments on one such model [^] lend 
support to this idea. In the inverse transport range the number of states decrease towards 
the energy dissipation range. A speculative conclusion is of this work is then that perhaps the 
true state of the 2D Navier-Stokes equation is a forward cascade of enstrophy, or a non-local 
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transfer with spectrum steeper than k~ 3 , on scales smaller than the forcing, but a formal 
statistical equilibrium on scales larger than the forcing. The latest numerical simulations by 
Borue [TjJ seem to agree with this picture. 
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FIGURES 

FIG. 1. The second order structure function u>2(n) = (j2j l^n,j| 2 y of shell vorticity as a func- 
tion of shell number, logarithmic scale. 

FIG. 2. The second order structure function of vorticity as a function of shell number in the 
range of forward transport of enstrophy, real scale. 

FIG. 3. Probability distribution function of the real parts of the shells variables W$ o an d W30 o 
in the range of forward transport of enstrophy. The full line is a Gaussian distribution with the 
same variance. 



FIG. 4. The kurtosis o>4(n)/o>2(n) 2 where 0*4(71) = ( \W n> j\ 2 ) ) of the shell vorticity dis- 



\ rr n,j 

tribution function. The horizontal line is the prediction if the three shell variables are assumed to 
be uncorrelated Gaussian variables with the same variance. 

FIG. 5. The mean enstrophy flux (triangle), and the mean squared enstrophy flux (square), in 
logarithmic scale. At the peak the standard deviation of the flux is about 50 times larger than the 
mean value. 

FIG. 6. The mean energy flux (triangle), and the mean squared energy flux (square), in loga- 
rithmic scale. At the peak the standard deviation of the flux is about 10, 000 times larger than the 
mean value. 

FIG. 7. Estimates of k% for a system of 50 shells from —24 to +25. The symbols refer to runs 
with forcing ( |l~6| ) v = 10~ 9 , v' = 10~ 6 and flux = 0.05 (square); v = 10~ 9 , v' = 10~ 6 and flux = 
0.1 (triangle); v = 10~ 6 , v' = 10 -4 and flux = 0.2 (full square). The full triangles refer to the data 
of figure ^. 

FIG. 8. Lyapunov exponent as a function of Reynolds number, defined as Re = u^ 1 . The 



system consists of 50 shells from —24 to +25 with forcing (16) with 7/ = 0.1 and v' = 10 . 
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FIG. 9. Spreading of a perturbation as a function of time. The bottom curve refers to time 
0.2 seconds from the perturbation. The others, bottom to top, after time intervals of 0.4 seconds. 
The system consists of 50 shells from —24 to +25 with forcing ( |i~6|) with rj = 0.1, v = 10 -9 and 
v' = 10 -6 . At the initial time the perturbation was of strength 10~ 10 . The curves are rather 
insensitive to the location of the initial perturbation, in figure it was located on shells —6 and —5. 
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